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Abstract 

In this paper, we examine a flat FLRW spacetime with a scalar field 
potential and show by applying Osgood’s criterion to the Einstein field 
equations that all such models, irrespective of the particular choice 
of potential develop finite-time singularities. That is, we show that 
solutions to the field equations become singular in finite time, which 
can have important implications for the role of inflation in cosmological 
models. We further point out that a possible reason for this behaviour 
is that the solutions to the field equations in such inflationary scenarios 
do not obey global existence and uniqueness properties, which is a 
typical characteristic of solutions that diverge in finite time. 


1 Introduction 

The purpose of this paper is to demonstrate the existence of future singular¬ 
ities that occur in finite time in inflationary models with an arbitrary scalar 
field potential within the context of a spatially flat Friedmann-Lemaitre- 
Robertson-Walker (FLRW) spacetime. The concept of the existence of a 
singularity in the context of general relativity has been studied quite ex¬ 
tensively. With respect to the cosmological case, there is the fundamental 
singularity theorem mmm which describes irrotational geodesic singular¬ 
ities, and states that if A < 0, p + 3p > 0, and p + p > 0 in a fluid flow 
where p is the energy density of the fluid, p is the pressure of the fluid, 
in addition to having u = 0, ui = 0, and Hq = 0 at some time so, then a 
spacetime singularity, where either the expansion scalar goes to zero or the 
shear scalar diverges occurs at a finite proper time To < I/Hq before so- 
Further elaborations are built upon this singularity theorem. It is of 
interest to note that in [1J five possible routes to avoid the conclusions of this 
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singularity theorem are discussed in detail. They are a positive cosmological 
constant, acceleration, vorticity, an energy condition violation, or alternative 
gravitational equations. We refer the interested reader to |lj for further 
discussions regarding these issues. 

Following [I], we also note that singularities occur in cosmology not only 
in the context of FLRW models, but also for realistic anisotropic and in¬ 
homogeneous models of the universe in which the strong energy condition 
p + 3p > 0 is satisfied. Related to this, Penrose [3j did pioneering work on 
black hole singularities producing a theorem that allowed one to predict the 
existence of singularities in realistic gravitational collapse cases. Hawking 
[5] extended these results to the cosmological context leading to the famous 
Hawking-Penrose singularity theorem which implied that space-time singu¬ 
larities are to be expected if either the universe is spatially closed or there 
is an ‘object’ undergoing relativistic gravitational collapse (existence of a 
trapped surface) or there is a point p whose past null cone encounters suf¬ 
ficient matter that the divergence of the null rays through p changes sign 
somewhere to the past of p (i.e. there is a minimum apparent solid angle, 
as viewed from p for small objects of given size). The theorem applies if 
the following four physical assumptions are made: (i) Einstein’s equations 
hold (with zero or negative cosmological constant), (ii) the energy density 
is nowhere less than minus each principal pressure nor less than minus the 
sum of the three principal pressures (the ‘energy condition’), (iii) there are 
no closed timelike curves, (iv) every timelike or null geodesic enters a region 
where the curvature is not specially aligned with the geodesic. (This last 
condition would hold in any sufficiently general physically realistic model.) 

Further, Bekenstein [6j studied exact solutions of Einstein-conformal 
scalar equations and presented a class of FLRW models which contained 
both incoherent radiation and a homogeneous conformal scalar field which 
bounce and never pass through a singular state, thereby circumventing the 
singularity theorems by violating the energy condition. Parker and Fulling 
pj considered a classical gravitational field minimally coupled to a quan¬ 
tized neutral scalar field possessing mass. They concluded that quantum 
effects can sometimes lead to the avoidance of the cosmological singularity, 
at least on the time scale of one Friedmann expansion. Collins and Ellis 
[8] examined in detail the singularities that occur in Bianchi cosmologies. 
Barrow and Matzner |9j showed that the singularity corresponding to ho¬ 
mogeneous and isotropic universes observationally equivalent to ours must 
be of simultaneous Robertson-Walker type containing only small curvature 
fluctuations. In an interesting application of Robertson-Walker singulari¬ 
ties, Barrow [TO] showed that the large-scale velocity and vorticity fields 
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required for the vortex and spinning-core theories of galaxy formation can 
be generated primordially in a natural way from a FLRW singularity. 

This paper in particular is concerned with the existence of singularities 
in inflationary models. Some work has been done in this regard. In particu¬ 
lar, Kohli and Haslam El proved the existence of finite-time singularities in 
the case of stochastic eternal inflation for an arbitrary inflaton field poten¬ 
tial. Nojiri, Odintsov, and Oikonomou |12] analyzed in detail the concept of 
singular inflation with scalar-tensor and modified gravity. Nojiri, Odintsov, 
Oikonomou and Saridakis m studied the existence of finite-time singu¬ 
larities in the inflationary and late-tinre regimes of cosmological evolution. 
Nojiri, Odintsov, and Tsujikawa m investigated in detail the properties 
of future singularities in a universe dominated by dark energy including 
a phantom-type fluid. Recently, Odintsov and Oikonomou m examined 
cases in which singular evolution can be consistently incorporated in defor¬ 
mations of the R 2 inflationary potential. Barrow and Graham m proved 
that a spatially homogeneous and isotropic universe containing a scalar field 
with potential V((f>) = Acj) n with 0 < n < 1 and A > 0 always develops a 
finite-time singularity at which the Hubble rate and its first derivative are 
finite, but its second derivative diverges. We are not aware of any other work 
in the literature that describes the existence of finite-time singularities in 
FLRW inflationary cosmologies. Note that throughout we used geometrized 
units where 8irG = c = 1. 


2 Description of the Model 


We consider a spatially homogeneous and isotropic universe described by 
the FLRW metric |1] as 

ds 2 = —dt 2 + a 2 (t) [dr 2 + f 2 (r) [d,0 2 + sin 2 ddq l 2 )] , u^ (1) 


where 

/(?’) = (sinh r, r, sinr) for K = (—1,0,+1), (2) 

where K denotes the sign of the curvature of the particular FLRW model 
under consideration. Namely, K = —1 refers to hyperbolic FLRW models, 
K = 0 refers to flat FLRW models, and K = +1 refers to positively curved 
FLRW models. 

We assume that energy-momentum tensor is that of a scalar field and 
has the form [1] 


Tf = V°0VV 


- V c 0V c (/» + v((j>) 


(3) 


3 




The Einstein field equations that also describe the dynamics of the model 
are given by the Raychaudhuri and Friedmann equations [I], 


H = -H 2 + - 
3 L 

12 


1 ( 3 ) 

3 H 2 = V{cj)) + ^- + - R, 


(4) 

(5) 


where H is the Hubble parameter, and ^ R is the three-dimensional Ricci 
scalar, which is constant for the FLRW models. 

The contracted Bianchi identities give an evolution equation for the 
scalar field, 0, which in this case is precisely the Klein-Gordon equation 

IB, 

4> + 3 Hj> + V'(<i>) = 0. (6) 

Together Eqs. 0, and 0 fully describe the dynamics of the cos¬ 
mological model described by Eqs. 0 and 0. 

We further define expansion-normalized variables denoted by X, Y as 


X = 




(7) 


As is done in mi, we additionally introduce a dimensionless time variable 
r, such that 


Substituting Eqs. 0 and 
dynamical system 


dt 1 
dr H 
into Eqs. 


( 8 ) 

), one obtains the following 


X' = X(l+q + YX), (9) 

y’ = -3X— ^(1 — Y 2 )+Y(l + q), (10) 


subject to the constraint 

X 2 + Y 2 = 1. (11) 

Note that q is the deceleration parameter found through Raychaudhuri’s 
equation and is in this case 


q = 2Y 2 — X 2 . 


( 12 ) 


Following m, we also have defined a model-dependent dimensionless 
parameter, 


A 


3 y 

2~V 


(13) 
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Note that in the definition of A in Eq. G3D. we have set the Planck mass 
to unity as per our choice of units. This will allow us to consider the dy¬ 
namics corresponding to an arbitrary inflaton potential making our results 
as general as possible. 

Making use of Eqs. (HID and ([HD in Eqs. @ and (fTUft . we see that the 
equations decouple, 

X' = X (3-3X 2 + , (14) 

Y' = (Y 2 — l) (3Y + A). (15) 

Notice that from Eq. m, that the energy-momentum tensor is that of 
a perfect fluid with energy density 

V- = + T(0), (16) 

and pressure 

p= l -^-V{ct>). (17) 

Following [I], we note that for a non-negative potential energy V > 0, one 
has the restriction that 

- < 1. (18) 

M 

For what follows, we will exclusively concentrate on the dynamics rep¬ 
resented by Eq. m- That is, we will analyze the initial value problem 

X' = X (3 - 3X 2 + y/l - a: 2 a) , 

X(0) = e (19) 

3 Osgood’s Criterion and Finite-Time Singulari¬ 

ties 

We will demonstrate the existence of finite-time singularities by making use 
of Osgood’s criterion applied to Eq. m Osgood’s criterion due to W.F. 
Osgood (1898) [T9] which is a very fundamental result describing conditions 
for the uniqueness of solutions to ordinary differential equations. As noted 
in m, Osgood established that the solution y(t ) of 

y '(t) = b (y(t)) , t> 0, (20) 

2 /( 0 ) = 

5 




diverges within a finite time, that is, y(t) —» ±oo if 

po ° ds 


s: 


< oo. 


( 21 ) 


/{ & ( s ) 

We further note that the time it takes for the solution to become singular 
is given by 


T = 


j: 


ds 


/ 5 b( s y 

such that 

T = sup {t > 0 : T < 00 } . 

Motivated by Osgood’s criterion, recalling Eq. (ED, we denote 


( 22 ) 

(23) 


F(X) = X (3 - 3X 2 + y/l - X 2 \) . (24) 

, and obtain 


We now compute the integral in (1211) using Eq. 

ro ° ds r°° dS 


r°° ds _ r 

.L ns) ' k 


5 5 (3 - 35 2 + Vl - 5 2 A) 

3 log 3 log 11 - V 1 - € 2 1 log 

A 2 - 9 2(3 + A) 


1 + V 1 - £ 2 3 log 3^/1 -e 2 + A 

■ + — 


2(3 - A) 


9-A 2 


(25) 


where X(0) = f. 

The integral as computed in Eq. (l25l) diverges if A = ±3. That is, per 
Osgood’s criterion and Eq. ED, there is no finite-time singularity if and 
only if 

(3 7'W 


2 v W) = ±3 ' (26) 

Except for the specific choice of the inflaton potential as described in 
Eq. (1261) . the types of inflationary models as described in this paper always 
admit a finite-time singularity, where X (r) —>• 00 , for r < 00 . 


4 Discussion 

We have shown thus far that inflationary models with arbitrary potential 
V(4>) in a flat FLRW background become singular in a finite time. Exploring 
this fact in more detail, recalling Eq. (1241) . we note that 


F'{X) =X 2 [- 


X 


Vi-x 2 


-9 + Xy/l - X 2 + 3. 


(27) 
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One therefore sees that F'(X ) is only continuous for any subinterval of 
— 1 < X < 1, albeit, it is also the domain for F'(X). This implies that 
F'(X) is also bounded over any closed sub interval of — 1 < X < 1, that 
is, F(X) is locally Lipschitz continuous over this domain. This implies by 
the existence and uniqueness theorem of ordinary differential equations m, 
that over any subinterval of — 1 < X < 1, solutions to Eq. (ED will exist 
and be unique. 

However, the computation of the integral in Osgood’s criterion in Eq. 
(|25l) also shows that while X(t ) —>• oo in finite time, it will first exit the 
domain of definition in finite time, that is, X(r) —> 1 in finite time. This 
is the point at which the singular solution begins to develop. Since, at this 
point, one can see that Eq. (1271) diverges, and so, all one has is the continuity 
of F(X) at X = 1, which only guarantees existence of solutions, but not 
uniqueness. For the Einstein field equations to be well-posed, one needs 
both existence and uniqueness of solutions, which is clearly not the case for 
the inflationary models described in this paper. 

5 Conclusions 

In this paper, we examined a flat FLRW spacetime with a scalar field poten¬ 
tial and showed through Osgood’s criterion that all such models, irrespective 
of the particular choice of potential develop finite-time singularities, which 
is due to the fact that solutions to the held equations rapidly diverge in a 
finite time. We further discussed that a possible reason for this behaviour is 
that the solutions to the held equations in inhationary models do not obey 
global existence and uniqueness, which is a typical characteristic of solutions 
that diverge in finite time. This, as it turns out, is a very important point. 

Even if solutions to a particular differential equation are unique over a 
hnite interval, one can continuously patch together such intervals to form 
a maximal interval on which the solution is dehned. That is, let ^ a (r) be 
the unique solution of the differential equation in question, which satishes 
V’a(O) = a i and let ( Tmin,T max ) denote the maximal interval on which ip a (t) 
is dehned. If t max is hnite, then [17] 

lim ||^ 0 (t)|| = Too. (28) 

T y Tmax 

That is, solutions can only be extended indehnitely unless the solution di¬ 
verges in a hnite time. However, we showed in this paper, that this is not 
the case, as solutions indeed diverge in a hnite time. 
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